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ABSTRACT 


An  approximate  analysis  is  performed  for  the  combustion  rate 
of  carbon  in  air  at  the  stagnation  region  of  a  blunt  body  in  hy¬ 
personic  flow  for  the  intermediate  regime  of  conditions  between 
the  reaction-kinetic-limited  and  the  diffusion- limited  extremes. 
The  effect  of  mass  addition  to  the  boundary  layer  caused  by  the 
combustion  is  considered.  Combustion  rates  in  the  intermedi¬ 
ate  regime  are  derived  using  closed-form  approximations  for 
available  boundary  layer  analyses.  The  results,  shown  in  a 
generalized  form,  are  valid  for  arbitrary  reaction  kinetics  and 
combustion  products.  The  limitations  of  the  series-resistance 
approach  are  demonstrated.  A  numerical  comparison  of  the 
closed- form  results  for  the  mass  loss  rates  in  the  diffusion- 
limited  regime  and  the  recent  results  of  Scala  indicate  that 
the  present  simplifications  are  satisfactory  for  engineering 
purposes. 
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I.  INTRODUCTION 


Ballistic  re-entry  vehicles  employing  charring-ablator  heat  shields  encounter 
high  surface  temperatures  during  re-entry  and  the  combustion  rate  of  the 
graphitic- char  material  is  limited  by  the  rate  at  which  oxygen  can  diffuse 
through  the  boundary  layer  to  the  burning  surface.  Lifting-types  of  re-entry 
vehicles,  however,  typically  encounter  lower  heat  transfer  rates  and  cor¬ 
respondingly  lower  surface  temperatures.  In  this  case,  a  graphitic -char 
materiad  or  a  graphite  may  combust  at  a  rate  that  is  more  or  less  dependent 
upon  the  heterogeneous  reaction-kinetics  of  the  surface  combustion  process. 
The  degree  to  which  the  combustion  is  dependent  upon  surface  kinetics  is  a 
function  of  the  flight  condition,  the  vehicle  configuration,  the  location  upon 
the  vehicle,  and  the  chemical  make-up  and  temperature  of  the  surface. 

In  1936,  Parker  and  Hottel  (Ref.  1)  reported  on  an  experimental  study  in 
which  the  combustion  rate  of  graphite  and  oxygen  was  measured  at  values  be¬ 
tween  the  reaction-kinetic  limit  and  the  diffusion  limit.  The  reaction-kinetic 
limit  may  be  defined  as  the  combustion  rate  under  conditions  where  the  oxy¬ 
gen  concentration  at  the  burning  surface  is  the  same  as  the  oxygen  concentra¬ 
tion  far  from  the  surface,  i.  e.  ,  at  the  edge  of  the  boundary  layer.  The  • 
diffusion- limit  is  defined  as  the  combustion  rate  under  conditions  where  the 
oxygen  concentration  at  the  burning  surface  is  essentially  zero  (for  practical 
ranges  of  wall  temperatures).  The  diffusion- limit  is  found  to  vary  only 
slightly  with  surface  temperature,  due  to  small  changes  in  boundary  layer 
transport  properties.  The  reaction-kinetic  limit  possesses  an  extreme  de¬ 
pendence  on  surface  temperature,  usually  expressed  as  an  exponential  after 
Arrhenius. 

Parker  and  Hottel  found  their  data  indicated  that  the  combustion  rate  of 
carbon  in  the  intermediate  range  of  surface  oxygen  concentrations  could  be 
expressed  by  a  relation  analogous  to  that  for  the  flow  of  electrical  current 


through  two  resistances  in  series.  They  identified  the  inverse  of  the 
reaction-kinetic  limited  combustion  rate  as  a  chemical  resistance  and  the 
inverse  of  the  diffusion- limited  combustion  rate  as  a  diffusion. resistance.  In 
1947,  Frank-Kamenetskii  (Ref.  2)  presented  a  simplified  analysis  of  combus¬ 
tion  in  the  intermediate  regime,  which  ignored  the  effect  on  the  boundary 
layer  characteristics  of  mass  transfer  from  the  wall.  According  to  Frank- 
Kamenetskii' s  analysis,  the  data  of  Parker  and  Hottel  conformed  to  fractional- 
order  rather  than  first-order  kinetics  as  they  had  originally  deduced. 

In  1961,  Moore  and  Zlotnick  (Ref.  3)  performed  an  approximate  analysis  of 
the  combustion  of  graphite  in  the  intermediate  regime  for  a  unity  Lewis  num¬ 
ber,  using  the  reaction-kinetic-data  of  Blyholder  and  Eyring  (Ref.  4).  They 
considered  the  case  of  a  boundary  layer  in  which  CO  evolving  from  the  heter¬ 
ogeneous  reaction  diffused  toward  a  gas  reaction  zone  at  some  arbitrary  dis¬ 
tance  from  the  surface  where  it  is  reacted  with  oxygen  to  form  CO^.  The 
CO^  then  diffused  back  to  the  graphite  surface  to  react  heterogeneously  with 
graphite  to  form  gaseous  CO. 

In  1962,  Scala  (Ref.  5)  reported  the  results  of  numerical  boundary  layer  cal¬ 
culations  for  the  combustion  rate  of  graphite  at  the  stagnation  point  of  a 
hemisphere  in  hypersonic  flow.  He  considered  the  cases  of  a  boundary  layer 
in  which  gas -phase  reactions  were  prohibited  (frozen)  and  a  boundary  layer 
in  which  gas -phase  reactions  proceeded  at  a  rate  sufficient  to  cause  local 
thermochemical  equilibrium  at  all  locations. 

In  the  intermediate  regime  (between  the  rate-limited  and  diffusion- limited 
combustion),  Scala  treated  the  kinetics  as  being  of  one-half  order  and  cal¬ 
culated  the  combustion  rate  by  the  series-resistance  method. 

Scala' s  numerical  auialysis  provides  an  accurate  solution  of  the  diffusion- limited 
combustion  rate  and  an  approximate  solution  of  the  intermediate -regime  com¬ 
bustion  rate.  His  intermediate  regime  results  must  be  considered  only  ap¬ 
proximate  since  it  has  not  been  shown  that  the  series-resistance  method 
utilized  may  be  applied  for  reactions  of  arbitrary  order,  different  from  unity. 
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In  this  report,  the  problem  analyzed  by  Scala  will  be  studied,  for  the  frozen 
boundary  layer  case,  by  the  use  of  approximate  but  closed-form  solutions  of 
the  stagnation- point  boundary  layer  equations.  The  purpose  of  the  analysis 
is  to  show  the  general  coupling  features  of  the  oxygen  diffusion,  the  surface  re¬ 
action  kinetics,  and  the  mass  transfer  of  combustion  products  with  regard  to 
the  overall  combustion  rate.  The  accuracy  of  the  series-resistance  method 
will  be  investigated.  In  addition,  it  will  be  shown  that  the  present  approxi¬ 
mate  results  can  predict  the  diffusion- limited  ablation  rates  within  satis¬ 
factory  accuracy. 
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II.  ANALYSIS 

The  problem  of  a  heterogeneous  reaction  at  a  wall  undergoing  a  transition 
from  a  dependence  upon  surface  kinetics  to  a  dependence  upon  boundary  layer 
diffusion,  for  the  case  of  wall- catalyzed  recombination  of  dissociated  atoms, 
was  treated  by  Chung  (Ref.  6)  in  1959.  In  that  analysis,  Chung  investigated 
the  effect  of  air  injection  through  a  porous  wall  of  arbitrary  catalytic  ef¬ 
ficiency  on  the  nature  of  the  transition  from  the  reaction-kinetic  limit  to  the 
diffusion  limit  of  wall -catalyzed  oxygen  recombination  rates.  A  direct  analogy 
between  the  graphite -combustion  problem  and  the  wall- catalyzed  atom  recom¬ 
bination  problem  exists  for  the  case  where  the  boundary  layer  conditions 
prohibit  gas-phase  reactions.  The  essential  features  of  the  ansdysis  of  Ref¬ 
erence  6  will  be  utilized  in  the  development  which  follows.  It  is  convenient 
to  choose  a.  flow  configuration  for  which  satisfactory  solutions  for  diffusion 
and  heat  transfer  are  available.  Such  a  configxiration  is  the  stagnation  region 
of  an  axisymmetric  or  two-dimensional  blunt  body  in  hypersonic  flow. 

The  fundamental  equations  for  a  frozen  laminar  boundary  layer  at  the  stagna¬ 
tion  region  sketched  in  Figure  1  are 

Continuity: 


(1) 


Momentum: 


+  pv 


/8u\  8  /  8uV 

\^)  =  '5y\^'Sy)  • 


8P 

■ST 


(2) 


Note:  c  a  0  for  a  two-dimensional  body,  1  for  a  rotationally- symmetric  body. 
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Figure  1.  Physical  Model 
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Energy: 
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'  Sy^-W)  ■  ■SyfVW  - 
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Diffusion: 

pu(^)tpv(|i)=|p^D|s) 

(4) 

The  following  transformation  relations  are  utilized: 

(  diji 

P“^o  =  -37 

(5a) 

c  ^ 

P'^’^O  = 

(5b) 

r*  Zf , 

•  =  X  PePe%"o 

(5c) 

II 

<n  ^ 

(5d) 

e 

(5e) 

It 

(5f) 

Where  C  is  assumed  to  be  constant  in  the  transformation.  Also, 

it  will  be 

assumed  that  Pr  and  Sc  numbers  are  constant  across  the  boundary  layer. 
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The  additional  quantities  g  and  Z  are  defined  as 


g(n)  =  ^ 


Z(ti)  = 


'l.e 


(5g) 


(5h) 


where  is  the  concentration  of  the  gaseous  reactant  of  interest,  oxygen. 
The  continuity  equation  is  satisfied  by  the  stream  function  transformation, 
and  the  momentum,  energy,  and  diffusion  equations  transform  to 


f.H  +  ffn  _  (f 

2 

=^(pjr- 


fZ’+^Z"  =0 


(6) 


(7) 


(8) 


Lees  (Ref.  7)  showed  that  the  terms  on  the  right-hand  sides  of  Eqs.  (6)  and 
(7)  can  be  neglected  without  appreciable  error  if  the  wall  surface  is  highly 
cooled.  This  condition  is  fulfilled  at  wall  temperatures  of  practical  interest 
for  graphite  materials  in  the  hypersonic  regime.  Thus  the  equations  simplify 
to 


f"  +  ff"  =  0 


fg'  +T^g"  =  0 


(9) 

(10) 
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The  boimdary  conditione  to  be  imposed  for  the  solution  of  the  momentum  and 
diffusion  equations  are 


f(0)  =  f^ 

(12a) 

f'(0)  =  0 

(12b) 

f'(oo)  =  1 

(12c) 

Z(oo)  =  1 

(12d) 

Z{0)  = 

(12e) 

It  is  noted  that  f^  and  are  not  arbitrary,  but  depend  on  the  surface 
chemical  reaction.  The  frozen  energy  equation,  (10),  will  not  be  solved,  but 
solutions  may  be  found  in  References  6  and  8. 

The  rate  of  mass  loss  from  a  burning  surface  is  expressed  as  the  Arrhenius 
form 


m  =  k, 
w  1 


■E  /RT 
a  w 


(13) 


where  k,  is  the  frequency  factor,  P,  is  the  partial  pressure  of  the  gaseous 
1  1  I  w 

reactant  at  the  wall,  and  n  is  the  order  of  reaction.  The  partial  pressure 
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P,  ^  can  be  expreeeed  in  terms  of  the  total  pressure  P  ,  molecular 

*  f  Wr 

weight  ratio  and  oxygen  reactant  concentration  at  the  wall  c^  ^ 


P 


l.w 


(14) 


The  frequency  factor  and  the  exponential  term  are  combined  as  k^. 

-E  /RT 

=  kj  e  ^  (15) 


Taking  the  static  pressure  as  constant  across  the  boundary  layer,  we  may 
express  the  mass  loss  rate  as 


m  s  k 
w  w 


w 


(16) 


The  stoichiometric  relation  for  heterogeneous  combustion  is  written  as 


n  M  +  n,  M,  -*•  n  M 
w  w  11  ^  p  p 


(17) 


where  w  refers  to  carbon,  1  to  oxygen,  and  p  to  the  product  species.  Thus 
the  gaseous  reactant  transport  rate  corresponding  to  the  wall  mass  loss  rate  is 


"'l  n  lA.  wl  e\M,  /  I  *^1, 
w  w  L  '  I  'wj 


w 


(18) 


Equating  this  to  the  net  transport  rate  of  gaseous  reactauit  at  the  wall,  &e 
boundary  condition  is  obtained  for  the  diffusion  equation  at  y  =  0 
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(19) 
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Using  the  relations 

ro{x)  =  X 

Pe^^e  =  ^PePe^ 

=  Px 

P^v^  =  +  ‘>^*PePe^s^ 

for  the  stagnation  region,  and  the  transformation  relations  (5),  the  diffusion 
boundary  condition  becomes 


Z'(0)  =  Sc 


“l^l  .  fp  /ia\  T  n-l 

tnsrKl^Jrssr)  ®i,e  .  i 

w  w  L  >  1  ^wJ  f  Z^'^fOv 


z”(0) 


(20) 


Since  the  wall  temperature  and  the  parameters  associated  with  it  are  assumed 

to  be  invariant  with  x,  similarity  is  not  impaired  by  the  inclusion  of  the 

heterogeneous  chemical  kinetics  in  the  boundary  condition  at  the  reacting 

surface.  It  is  convenient  to  define  a  DamkShler  number  D  as 

m 


D 


m 


n,M, 

n  ^  w 
w  w 


>/rr+TTPip^{r^Ts^ 


1 


(21) 
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This  parameter  characterise*  the  ratio,  for  a  unit  maaa  of  oxygen,  of  the 
diffuaion  time  at  the  diffuaion  limit  to  the  reaction  time  at  the  reaction- 
kinetic  limit,  and  allow*  the  following,  more  compact,  veraion  of  equation  (20): 

Z'(0)  =  Sc[d^  -  f^  Z^-”(0)]z"{0)  (22) 

The  application  of  boundary  condition  (12d)  to  the  diffuaion  equation  (11) 
lead*  to  the  aolution 


1  -  Z(0)  = 


dq 


(23) 


For  low  value*  of  f^,  the  following  relations  furnish  good  approximations  of 
the  numerical  results  of  References  6  and  9  for  air  boundary  layers: 


P” 


(T1)]®®  dTiaC, 


f  <  0.  2 
w 


(24) 


[f"(0)]^‘=arC2+  €3!^ 


f  <  0.  2 

w 


(25) 


where  for  Sc  =  0.  72  (Ref.  8)^ 

Cj  =  1. 391 


C2  =  0.  580 


and  for  Sc  =  0.  514  (Ref.  6) 


Cj  =  1.  825 


C2  =  0.  678 


C3  =  0.610 


C3  =  0.485 


These  are  Low's  results,  numerical^  modified  to  correspond  to  the  present 
nomenclature.  There  is  a  factor  of  '>/2  involved  between  the  present  f  and  q 
and  those  of  Low.  These  factors  are  pointed  out  by  Chung  (Ref.  6). 
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When  an  injected  epecie  has  a  molecular  weight  much  different  from  that  of 
air,  these  results  must  be  modified:  however,  the  CO  and  CO2  combustion 
products  for  this  problem  have  molecular  weights  sufficiently  close  to  air  that 
such  modifications  would  be  small.  Also,  the  f^  values  in  this  problem  are 
characteristically  below  0.  2  even  at  the  diffusion  limit.  Therefore,  these 
approximations  are  considered  applicable  within  reasonable  accuracy. 

By  applying  the  boundary  condition  at  the  wall  (Eq.  22)  to  the  diffusion  equation 
solution  (Eq.  23)  and  expressing  f^  in  terms  of  the  Damkdhler  parameter, 
we  obtain,  after  rearranging  terms. 


Now,  for  the  given  values  of  Sc  and  D^,  and  for  the  given  products  of  com¬ 
bustion,  Eq.  (26)  gives  the  wall  concentration  of  oxygen.  The  combustion 
rate  can  then  be  obtained  from  Eqs.  (16)  and  (21). 

The  diffusion- limited  combustion  rate  may  be  calculated  by  letting  Z(0)  0  in 

Eq.  (26)  and  solving  for  D^Z^(O).  Then  the  combustion  rate  is 


n  M 
w  w 


'  w,d.  1. 


(27a) 


or 


(27b) 


lU.  RESULTS  AND  COMPARISON  WITH  OTHER  INVESTIGATORS 


The  first  point  to  be  examined  ie  that  of  the  accuracy  of  this  method  of 
calculating  the  diffusion-limited  mass  loss  rate.  The  recent  results  of 
Scala  (Ref.  5),  obtained  using  accurate  numerical  solutions  of  the  boundary 
layer  differential  equations,  afford  a  good  reference.  It  will  be  assumed  that 
the  ratio  of  carbon  monoxide -to -carbon  dioxide  gas  produced  by  the  combus¬ 
tion  is  given  by  the  equilibrium  relation. 

C(a)  +  C02(g)  2CO(g) 

The  equilibrium  constant  for  this  relation,  as  used  by  Scala,  is 


•Pco'" 


.  8.  80 
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Lees  (Ref.  7)  has  suggested  that  the  Schmidt  number  for  boundary  layers 
adjacent  to  burning  carbon  surfaces  is  between  0.  5  and  0.  7.  Therefore,  the 
relations  approximating  the  results  of  References  6  and  9  for  Schmidt  numbers 
of  0.  514  and  0.  72,  respectively,  furnish  useful  solutions  for  comparison. 


Figure  2  shows  the  diffusion-limited  mass  loss  rate  resulting  from  the 
solution  of  Eq.  (27b)  for  the  following  flight  conditions  of  Reference  5: 
vehicle  velocity,  20, 000  ft /sec;  altitude,  100,000  ft.  The  numerical  results 
of  Reference  5  are  also  shown  in  Figure  2.  A  large  increase  in  mass  loss 
rate  occurs  as  wall  temperature  is  raised  from  1000*^R  to  2000*^R,  due  to 
the  change  of  the  equilibrium  product  from  primarily  CO^  to  primarily  CO 
gas.  The  small  variation  in  mass  loss  i^ate  occurring  at  temperatures 
higher  than  2000°R  is  due  to  changes  in  the  density  and  viscosity  in  the 
boundary  layer. 
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The  results  of  Scale  correspond  to  the  solution  for  Sc  3  0.  514  at  low  wall 
temperatures  and  to  a  solution  that  would  result  using  a  Schmidt  number 
midway  between  0.  514  and  0.  72  at  the  higher  wall  temperatures.  It  is 
apparent  that  the  simplified  closed>form  solutions  can  yield  satisfactory 
engineering  results  for  the  diffusion -limited  mass  loss  rates  since  errors 
up  to  fifteen  percent  (Ref.  5),  due  to  uncertainties  in  transport  properties, 
are  possible  in  all  results  shown  in  Figure  2. 

In  the  intermediate  regime,  it  is  not  necessary  to  specify  particular  numeri¬ 
cal  values  for  the  Arrhenius  relation  since  the  Oamkohler  number  includes 
this  factor  in  the  general  results.  It  is,  however,  necessary  to  specify  the 
order  of  reaction,  n,  and  the  product  of  the  reaction.  As  discussed  by 
Scala  (Ref.  5),  there  is  considerable  disagreement  in  the  literature  on  com¬ 
bustion  as  to  the  order  of  reaction  for  -graphite  combustion.  One-half 
order  kinetics  are  widely  utilized,  however,  and  the  Scala  (Ref.  5)  and 
Moore  and  Zlotnick  (Ref.  3)  articles  both  incorporated  that  order.  Hottel 
found  experimentally  that  the  order  was  close  to  unity.  The  effective  com¬ 
bustion  product  can  range  from  primarily  CO  to  primarily  CO^.  To  accom¬ 
modate  this  range  of  reaction  orders,  products,  and  the  expected  variation 
of  Schmidt  number,  £kj.  (26)  was  solved  for  n  s  1  and  1/2;  CO  and  CO^ 
products:  and  Sc  s  0.  514  and  0.  72. 

The  results  of  this  calculation,  shown  in  Figure  3,  indicate  the  dimension¬ 
less  concentration  of  oxygen  at  the  wall  versus  the  Damkohler  number.  The 
average  molecular  weight  of  the  gas  mixture  at  the  reacting  surface,  (^)^> 
was  taken  as  that  of  air,  since  nitrogen  is  the  predominant  specie  and  the 
molecular  weights  of  the  combustion  products  are  not  far  different  from 
oxygen.  The  resume  for  a  CO^  product  are  shown  only  as  perturbations 
from  the  CO  product  curves.  The  differences  between  CO2  and  CO  product 
results  are  small  for  given  reaction  orders  and  Schmidt  numbers.  The 
present  results  for  unity  reaction  order  are  seen  to  require  an  increase  of 
four  orders  of  magnitude  in  IhunkShler  number  in  order  to  lower  the  oxygen 
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Figure  3.  Oxygen  Concentration  at  the  Wall  in  the  Intermediate  Regime 


concentrations  at  the  wall  from  99  percent  to  1  percent  of  the  free-stream 
concentration.  A  reaction  order  of  one -half  requires  ati  increase  in  of 
only  three  orders  of  magnitude  for  a  similar  transition.  The  different 
Schmidt  numbers  at  a  given  reaction  order  show  results  converging  at  the 
extremes  and  differing  a  small  amount  at  intermediate  DamkShler  numbers. 
Results  calculated  from  relations  derived  in  the  Appendix  using  the  series - 
resistance  method  are  also  shown  in  Figure  3  at  particular  DamkShler 
values.  It  is  noted  that  the  series -resistance  method  leads  to  oxygen  con¬ 
centrations  that  are  approximately  correct  for  unity  reaction  order  but 
quite  inaccurate  for  a  one-half  reaction  order.  The  series  results  are  not 
exactly  valid  even  for  unity  reaction  order  because  the  modifications  to  the 
boundary  layer  diffiision rates  due  to  mass  addition  are  not  adequately 
represented  in  the  intermediate  regime. 

The  fractional  mass  loss  rates,  defined  as  the  ratio  of  the  wall  mass  loss 
rate  to  the  wall  mass  loss  rate  at  the  diffusion  limit,  were  obtained  from 
the  results  of  Figure  3  and  Eqs.  (16)  and  (27b),  and  are  shown  in  Figure  4 
for  a  CO  product.  The  mass  loss  rates  given  by  the  reaction-kinetic  limit  and 
bythe  series  method  are  shown  for  reference.  The  present  results  require 
changes  of  the  same  order  as  in  Figure  3  to  accommodate  a  l-to-99 
percent  transition,  and  again  the  series  result  is  adequate  for  n  =  1  and 
inadequate  for  n  =  1/2  reaction  order. 

The  fundamental  reason  for  the  failure  of  the  series -resistance  method  when 
utilized  for  a  reaction  order  of  one -half  is  that  the  reaction-kinetic  limit 
then  does  not  depend  linearly  upon  the  oxygen  concentration,  whereas  the 
diffusion -limit  does  show  a  linear  dependence.  A  series -resistance  calcula¬ 
tion  is  predicated  upon  the  quantity  of  flow  having  a  linear  dependence  on  the 
potential  across  each  resistance. 

There  exist  considerable  uncertainties  in  the  values  of  the  reaction-kinetic - 
constants  for  the  combustion  of  the  various  grades  of  graphites  with  oxygen. 


Figure  4.  Fractional  Mass  IjOSS  Rate  in  the  Intermediate  Regime 


and  in  the  proportions  of  CO  and  CO^  produced  by  the  reaction.  It  is  felt 
that  the  present  results  are  sufficiently  simple  and  general  to  furnish  a 
useful  method  for  calculating  combustion  rates  with  any  kinetics  and  products 
that  may  be  found  appropriate.  The  question  of  the  proportions  of  CO  and 
CO2  produced  in  the  intermediate  regime  is  at  present  unanswered.  Although 
the  equilibrium  written  between  carbon,  carbon  dioxide,  and  carbon  monoxide 
is  a  reasonable  postulate  for  the  diffusion-limited  condition,  it  may  become 
invalid  in  the  intermediate  regime.  The  reaction-kinetics  of  the  carbon- 
carbon  dioxide  reaction  (as  described,  for  instance,  in  Ref.  2)  may  then 
become  significant.  Thus,  the  numerical  calculation  of  combustion  rates 
in  the  intermediate  regime,  at  the  present  time,  must  beperformed  for  the 
limiting  cases  of  either  CO  or  CO2  product  formation  to  show  the  possible 
ranges  of  these  results. 
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IV.  CONCLUSION 


The  simultaneous  effects  of  heterogeneous  reaction -kinetics  and  oxygen 
diffusion  to  a  burning  graphite  stagnation  surface  may  be  described  by  a 
simple  relation  (Eq.  26),  which  utilizes  closed-form  approximations  to 
available  boundary  layer  diffusion  solutions  and  a  DamkShler  number.  The 
OamkShler  number  characterizes  the  ratio  of  the  time  for  a  quantity  of 
oxygen  to  diffuse  to  the  surface  to  the  time  for  a  quantity  of  oxygen  to  chem¬ 
ically  react  with  the  surface. 

The  use  of  closed-form  approximations  to  available  boundary  layer  solutions 
allows  rapid  estimation  of  diffusion-limited  combustion  rates,  within  reason¬ 
able  engineering  accuracy. 

The  relations  derived  for  the  oxygen  concentration  at  the  burning  graphite 
surface  in  the  intermediate  regime  (Eq.  26)  may  be  utilized  for  any  reaction 
order,  reaction-kinetic  constants,  or  reaction  products  that  may  be  consid¬ 
ered  appropriate. 

The  proportions  of  carbon  monoxide  and  carbon  dioxide  products  formed  by 
combustion  in  the  intermediate  regime  are  presently  in  question;  therefore, 
the  numerical  value  of  the  carbon  mass  loss  rate  has  been  calculated  only 
for  the  limiting  cases  of  CO  or  CO^  product  formation.  An  investigation 
is  now  in  progress  to  extend  this  analysis  to  determine  the  effective  ratio 
of  carbon  monoxide  to  carbon  dioxide  produced  by  the  heterogeneous  reaction. 

A  series -resistance  approach  to  the  intermediate  regime  is  adequate  for  unity 
reaction  order  providing  the  diffusion -limited  mass  loss  rate  takes  proper 
account  of  the  effects  of  mass  addition.  The  series  approach  is  inadequate 
for  a  reaction  order  of  one -half. 
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where  r.  1.  refers  to  the  reaction-kinetic  limit  where  Z(0)  s  1  and  d.  1.  refers 
to  the  diffusion  limit  where  Z(0)  is  essentially  zero.  When  the  relations 
derived  in  Eqs.  (16)  and  (27b)  are  inserted,  these  terms  become 


Rearranging, 


1 


By  definition  of  th«  Demkflhler  number  (Eq.  2 1),  this  becomes 


or 
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This  derivation  has  been  made  only  for  the  purpose  of  allowing  a  comparison 
with  the  present  results.  As  shown  in  the  main  body  of  this  report,  the 
series  method  leads  to  erroneous  results. 
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NOMENCLATURE 


c 


C 


^r^2*^3 

c 

p 

o 


g 

h 

k 


k 

m, 

M 


n 

P 

Pr 

r 

o 

R 

■ 


maas  fraction  concentration 

Constanta  defined  in  Eqs.  (24)  and  (25) 
specific  heat 
diffusion  coefficient 

modified  OamkShler  number  defined  in  Eq.  (21) 
activation  energy 

transformation  variable,  -p^v^/y/il  +  c  )P(p^|i^  )^C 

dimensionless  enthalpy 

frozen  enthalpy 

thermal  conductivity 

frequency  factor 

specific  reaction  rate 

mass  loss  rate  of  surface 

molecular  weight 

reaction  order  or  stoichiometric  coefficient 

static  pressure,  atm 

Prandtl  number 

radial  coordinate 

universal  gas  constant 

body  nose  radius 

transformation  variable 
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NOMENCLATURE  (Continued) 


Sc 

Schmidt  number,  p/pD^^ 

T 

temperature,  absolute 

u 

velocity  component  parallel  to  wall 

V 

velocity  component  normal  to  wall 

X 

coordinate  parallel  to  wall  surface 

y 

coordinate  normal  to  wall  surface 

z 

oxygen  concentration  ratio 

P 

external  free  stream  velocity  gradient 

( 

exponent  for  body  shape 

n 

similarity  variable 

viscosity 

p 

density 

'I' 

stream  function 

Subscripts  and  Superscripts 

e 

edge  of  boundary  layer 

w 

at  the  wall  surface 

D 

average 

1 

gaseous  reactant 

d.l. 

diffusion  limit 

r.  1. 

reaction-kinetic  limit 

00 

vehicle 
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